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Abstract 

In  the  concept  known  as  Magnetized  Target  Fusion 
(MTF)  in  the  United  States  and  Magnitnoye  Obzhatiye 
(MAGO)  in  Russia,  a  preheated  and  magnetized  target 
plasma  is  hydrodynamically  compressed  to  fusion 
conditions.  Because  the  magnetic  field  suppresses  losses 
by  electron  thermal  conduction  in  the  fuel  during  the 
target  implosion  heating  process,  the  implosion  velocity 
may  be  much  smaller  than  in  traditional  inertial 
confinement  fusion.  Hence  “liner-on-plasma” 
compressions,  magnetically  driven  using  relatively 
inexpensive  electrical  pulsed  power,  may  be  practical. 
The  relatively  dense,  hot  target  plasma,  with  starting 
conditions  O(1018  cm'3,  100  eV,  100  kG),  may  spend  10 
or  more  microseconds  in  contact  with  a  metal  wall  during 
formation  and  compression.  Influx  of  a  significant 
amount  of  high-Z  wall  material  during  this  time  could 
lead  to  excessive  cooling  by  dilution  and  radiation  that 
would  prevent  the  desired  near- adiabatic  compression 
heating  of  the  plasma  to  fusion  conditions. 
Magnetohydrodynamic  (MHD)  calculations  including 
detailed  effects  of  radiation,  heat  conduction,  and 
resistive  field  diffusion  are  being  done,  using  several 
different  computer  codes,  to  investigate  such  plasma-wall 
interaction  issues  in  ongoing  MTF  target  plasma 
experiments  and  in  proposed  liner-on-plasma  MTF 
experiments. 

I.  INTRODUCTION 

Magnetized  Target  Fusion  (MTF)  is  an  approach  to 
controlled  fusion  that  is  intermediate  between  magnetic 
confinement  and  inertial  confinement  fusion  (ICF)  in 
time  and  density  scales.  Bigger  targets  and  much  lower 
initial  target  densities  than  in  ICF  can  be  used,  reducing 
radiative  energy  losses.  Reduced  losses  permit  near- 
adiabatic  compression  of  the  fuel  to  ignition 
temperatures,  even  at  low  (e.g.,  1  cm/psec)  implosion 
velocities.  In  MTF,  the  convergence  ratio  (rtaiUal  /rnMl)  of 
the  pusher  in  quasi-spherical  geometries  may  potentially 
be  less  than  10,  depending  upon  the  initial  temperature  of 
the  fuel  and  the  adiabaticity  of  the  implosion. 
Therefore,  “liner-on-plasma”  compressions,  magnetically 
driven  using  relatively  inexpensive  electrical  pulsed 
power,  may  be  practical  [1-4]. 

An  MTF  system  requires  two  elements:  (1)  a 
preheated  and  magnetized  initial  “target”  plasma;  (2)  a 
target  implosion  driver.  Because  the  reduced  energy 


losses  in  MTF  relax  the  power  and  intensity  requirements 
for  an  implosion  driver,  an  optimal  driver  source  for 
MTF  might  be  relatively  inexpensive  electrical  pulsed 
power,  to  drive  a  liner-on-plasma  implosion.  This  could 
utilize  either  fixed  pulsed-power  facilities,  such  as  the 
Pegasus  and  Atlas  capacitor  banks  at  Los  Alamos  and 
Shiva-Star  capacitor  bank  at  the  Air  Force  Weapons 
Laboratory  (Albuquerque),  or  explosive-flux- 
compression  generators,  such  as  Los  Alamos’  Procyon  or 
the  Russian  200-MJ-class  disk  flux  compression 
generators  [3,4],  Such  energy-rich  sources  might  allow  a 
demonstration  of  fusion  ignition  via  MTF,  without  a 
major  capital  investment  in  driver  technology. 

The  success  of  magnetized  target  fusion  will  depend 
upon  a  number  of  issues.  The  initial  target  plasma  must 
meet  minimum  temperature  (-50  eV,  preferably  100-300 
eV),  density  (between  10'3  andlO'6  g/cm3),  and  magnetic 
field  (>50  kG)  requirements,  and  must  have  a  lifetime, 
adjacent  to  the  supporting  wall,  greater  than  the 
implosion  time  (typically  several  psec  for  a  pulsed- 
power-driven  implosion).  Plasma-wall  interaction  must 
not  create  dynamical  effects  or  excessive  introduction  of 
impurities,  which  might  lead  to  rapid  cooling  of  the 
plasma.  The  target  plasma  must  be  readily  integrable 
with  drivers  for  compression  to  fusion  conditions.  As  it 
implodes,  the  liner  must  remain  sufficiently  intact  that  it 
can  effectively  compress  the  target  plasma. 

Los  Alamos  National  Laboratory  (LANL)  and  the 
All-Russian  Scientific  Research  Institute  of  Experimental 
Physics  (VNIIEF)  have  pursued  MAGO/MTF  research  in 
recent  years  [3-12].  Los  Alamos  is  presently 
investigating  three  candidate  target  plasmas:  the  Russian- 
originated,  explosive-pulsed-power-driven  MAGO 
plasma  formation  scheme,  the  high-density  Z-pinch,  and 
the  Field  Reversed  Configuration  (FRC),  an  elongated 
compact  toroid.  MAGO  work  includes  ongoing  joint  US- 
Russian  experiment  and  theory  aimed  at  determining  the 
suitability  of  the  plasma  created  for  MTF  compression. 
A  partially  wall-supported  deuterium-fiber-initiated  Z- 
pinch  experiment  at  LANL  has  been  investigated  for 
MTF  applications.  Los  Alamos  is  now  beginning 
experimental  and  theoretical  investigation  of  an  FRC 
plasma  for  MTF  compression. 

Related  experimental  and  computational  work  aims 
to  evaluate  explosive-flux-compression  generators  and 
existing  pulsed  power  facilities  as  MTF  liner  drivers.  A 
joint  LANL-VNIIEF  experiment  (high  energy  liner 
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“HEL-1”)  executed  in  August,  1996,  put  approximately 
100  MA  through  a  24-cm-radius,  4-mm-thick,  10-cm- 
long  aluminum  cylindrical  liner,  which,  by  the  time  it 
reached  the  5.5  cm  radius  of  the  central  measuring 
(diagnostic)  unit  (CMU),  was  imploding  at  7  mm/psec 
and  contained  approximately  20  MJ  of  kinetic  energy 
[7,8,11,12].  Such  a  liner  approaches  the  energy  and 
velocity  regime  required  to  drive  an  MTF  target  such  as 
MAGO  to  fusion  conditions.  Recent  liner  experiments  at 
Shiva-Star  (e.g.,  Faehl,  et  al.,  reported  at  this  conference) 
also  demonstrate  liner  performance  suitable  for  some 
MTF  target  plasmas. 

In  this  paper,  we  describe  computational  modeling  of 
proposed  cylindrical  liner-on-plasma  compressions  of 
such  target  plasmas  with  such  liner  systems. 
Computational  models  of  the  MAGO  target  plasma  and 
the  HEL-1  liner  have  shown  good  agreement  with 
experiment  [5-12],  These  computational  models  provide 
good  starting  points  for  liner-on-plasma  simulations.  The 
results  of  such  liner-on-plasma  simulations  can  point  out 
important  issues  which  must  be  confronted  in  designing 
such  experiments. 

II.  COMPUTATIONS 

Two  Los  Alamos  MHD  codes  have  been  employed 
in  modeling  MTF  experiments.  The  MHRDR  code  has 
been  used  extensively  in  modeling  the  MAGO  and  fiber- 
initiated  Z-pinch  target  plasma  experiments  (13). 
Because  MHRDR  is  fully  implicit  and  has  a  generalized 
Eulerian  structure,  in  which  zonal  quantities  such  as 
velocity  are  computed  relative  to  a  pre-programmable 
grid  velocity,  the  code  can  run  liner-on-plasma 
calculations  relatively  quickly.  However,  it  can  only 
compute  a  single  material  (e.g.,  the  target  plasma  DT); 
the  outer  radial  boundary  of  the  computational  mesh  is 
programmed  to  implode  at  the  expected  liner  velocity. 
Boundary  conditions  are  idealized:  perfectly  electrically 
conducting  and  zero-temperature,  infinitely  heat- 
conducting  walls  (hence  the  heat  conductivity  of  the 
magnetized  plasma  is  the  only  impediment  to  heat  flow 
out  the  boundary).  In  these  cases  MHRDR  uses 
cylindrical  r-z  geometry  with  a  single  Be  magnetic  field 
component.  The  plasma  radiative  energy  losses  and 
electrical  resistivity  are  taken  from  Los  Alamos 
“Sesame”  data  tables.  Heat  conduction  is  full  (arbitrary 
coota)  Braginskii,  unless  we  override  this  for  comparison 
to  other  models. 

A  second  Los  Alamos  Eulerian  MHD  code  [14]  was 
used  to  model  the  HEL-1  and  Shiva-Star  liner 
experiments.  This  code  can  compute  multiple  materials, 
such  as  an  aluminum  liner  imploding  onto  a  DT  plasma. 
However,  it  runs  much  slower  than  MHRDR  because  it  is 
not  fully  implicit  and  must  use  a  fixed  grid;  to  date  it  has 
only  been  possible  to  do  one-dimensional  liner-on-plasma 
calculations  (two-dimensional  runs  are  planned).  At 
present,  Braginskii  heat  conduction  is  only  included  for 
the  electron  fluid  (ion  fluid  has  the  non-magnetized 
conductivity  value);  however,  it  has  been  possible  to 


estimate  the  Braginskii  effect  on  the  ions  and  include  a 
constant  factor  times  the  non-magnetized  value  (1/15)  to 
approximate  the  full  Braginskii  magnetoinsulation  effect. 
Braginskii  ion  heat  conduction  is  being  added  to  this 
code.  The  diffusion  of  magnetic  field,  heat,  and  radiation 
between  liner  and  plasma  is  computed,  but  the  boundary 
between  liner  and  plasma  remains  sharp  (i.e.,  no 
intermixing  takes  place). 

A  one-dimensional  liner-on-plasma  problem  based 
on  demonstrated  target  plasma  and  liner  drive  quantities 
was  run  with  both  codes.  The  target  plasma  resembled  a 
late-time  (smaller  chamber)  MAGO  plasma:  100  eV,  10'5 
g/cm3,  with  3.0  MA  on  the  1.0  cm-radius  copper  inner 
conductor,  with  the  outer  aluminum  liner  starting  at  5.4 
cm  inner  radius  (0.6  cm  thick).  The  current  driving  this 
liner  was  based  on  a  portion  of  the  measured  HEL-1 
current,  which  would  implode  this  liner  from  5.4  cm  to 
1.36  cm  in  7.5  psec,  with  a  final  implosion  velocity  of 
1.6  cm/psec  (computed  liner  inner  radius  vs.  time  values 
were  used  in  the  MHRDR  calculation).  Such  a  liner 
implosion  is  not  an  optimized  choice  for  an  experiment 
intended  to  achieve  fusion  conditions,  but  it  represents 
something  clearly  achievable,  and  which  could  serve  as  a 
useful  step  in  demonstrating  compression  heating  of  an 
MTF  plasma. 

For  an  equivalent  volumetric  adiabatic  compression 
(cylindrical),  the  temperature  of  the  gas  would  go  from 
100  eV  to  1.03  keV.  The  MHRDR  one-dimensional, 
two-temperature  calculation,  which  includes  radiative 
and  conductive  heat  losses,  reaches  a  mass-weighted 
average  temperature  (<Ti>~<Tc>)  of  620  eV,  with  a  final 
radial  temperature  profile  peaking  at  900  eV  (Figure  1). 
The  final  radial  density  profile  is  relatively  flat,  about 
3X104  g/cm3,  except  for  higher  values  adjacent  to  the 
cold  walls,  particularly  the  imploding  wall.  Pressure  is 
peaked  similarly  to  temperature,  while  magnetic  field  is 
the  inverse  of  this.  A  MHRDR  three- temperature 
calculation,  starting  with  a  radiation  temperature  of  1  eV, 
with  open  boundaries,  duplicates  this  result,  confirming 
our  assertion  that  the  DT  plasma  is  optically  thin,  so  that 
two-temperature  calculations  with  a  radiative  energy  loss 
term  are  sufficient.  A  two-temperature  MHRDR 
calculation,  with  Braginskii  electron  heat  conduction  and 
non-magnetized  ion  heat  conduction  multiplied  times  an 
arbitrary  factor  of  1/15,  also  gave  a  final  peak 
temperature  of  900  eV;  this  multiplier  was  then  used  in 
calculations  with  the  other  code  to  approximate  the 
Braginskii  ion  insulation  effect. 

Calculations  with  the  second  code  show  similar  final 
plasma  conditions  in  the  target  plasma,  although  the  peak 
temperature  reached  is  slightly  under  800  eV,  compared 
to  the  idealized  MHRDR  case  900  eV.  Figures  2  and  3 
show  radial  temperature  and  current  distribution  in  the 
inner  rod,  target  plasma,  and  liner,  at  a  time  slightly  later 
than  the  MHRDR  profiles  in  Figure  1.  Note  the 
substantial  diffusion  of  heat  and  field  into  the  inner  and 
outer  metal  walls,  leading  to  wall  temperatures  as  high  as 
40  eV.  Also  note  that  plasma  pressure  is  high  enough 
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that  the  inner  rod  has  been  compressed  from  its  initial 
radius  of  1.0  cm  to  about  0.93  cm. 

MHRDR  has  also  been  used  to  compute  one-  and 
two-dimensional  compressions  of  computed  late-time 
MAGO  target  plasmas  at  2  cm/psec.  The  starting 
conditions  for  these  calculations  were  computed  second- 
chamber  plasmas  at  12  psec  in  the  LANL-VNIIEF 
MAGO-2  target  plasma  experiment;  the  MAGO-2  target 
plasma  calculations  showed  good  agreement  to  the  earlier 
time  experimental  measurements  available  [7-11].  The 
computed  late-time  MAGO  plasma  profiles  are  diffuse, 
wall-supported  Z-pinch  equilibria  which  show 
Kadomtsev  stability  to  m=0  perturbations.  In  liner-on- 
plasma  calculations,  the  10-cm  outer  wall  was  imploded 
at  2  cm/psec  to  a  final  radius  of  1.4  cm  (inner  wall  was 
1.2  cm).  One-dimensional  calculations  reached  a  peak 
mass-weighted  average  temperature  of  4.75  keV,  with  a 
peak  profile  temperature  of  7  keV.  Two-dimensional 
calculations  have  been  run  as  far  as  3.5  psec  to  date,  with 
the  <T>  the  same  as  in  the  one-dimensional  result  (800 
eV).  An  interesting  feature  can  be  seen  in  the  two- 
dimensional  calculation:  formation  of  convective  cooling 
cells  close  to  the  imploding  outer  boundary  (Figure  4). 
The  MHRDR  code  has  previously  been  used  in  studies  of 
such  cells  [15-17].  Since  the  bulk  temperature  reached 
has  not  changed  compared  to  the  one-dimensional  case, 
this  appears  (to  the  time  calculated  to  date)  to  be  a 
localized  effect  countered  by  stronger  heating  and 
insulating  processes. 

Will  ionized  wall  material  mix  with  and  cool  DT 
plasma  before  it  can  be  compressively  heated  to  fusion 
conditions  and  produce  significant  fusion  energy?  The 
answer  to  this  question  depends  upon  the  rates  of  the 
competing  processes  of  implosion,  heating,  mixing,  and 
cooling.  Detailed  calculations,  which  must  ultimately  be 
validated  by  experiment,  can  answer  this  question.  The 
codes  described  here  contain  substantial  portions  of  the 
physics  governing  these  competing  processes,  including 
potential  two-dimensional  effects.  Guided  by 
experimental  data  as  it  becomes  available,  we  can  utilize 
these  tools  to  predict  the  important  issues  for  optimizing 
the  design  of  future  MTF  liner-on-plasma  demonstration 
experiments. 

III.  CONCLUSIONS 

Magnetized  Target  Fusion  (MTF)  is  an  approach  to 
controlled  fusion  which  potentially  avoids  the  difficulties 
of  the  traditional  magnetic  and  inertial  confinement 
appoaches.  It  appears  possible  to  investigate  the  critical 
issues  for  MTF  at  low  cost,  relative  to  traditional  fusion 
programs,  utilizing  pulsed  power  drivers  much  less 
expensive  than  ICF  drivers,  and  plasma  configurations 
much  less  expensive  than  those  needed  for  full  magnetic 
confinement.  Computational  modeling  of  separate  MTF 
target  plasma  and  liner  implosion  experiments  has  shown 
good  agreement  to  experiment.  Combining  target  plasma 


and  liner  implosion  computational  models  allows  detailed 
theoretical  investigation  of  important  issues  for  proposed 
MTF  liner-on-plasma  experiments. 
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Figure  1.  Temperature  (eV)  vs.  radius  (m.),  MHRDR  1-d 
calculation. 
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Figure  2.  Temperature  (eV)  vs.  radius  (cm.),  second  code 
1-d  calculation. 


Figure  3.  Current  (MA)  vs.  radius  (cm.),  second  code  1-d 
calculation. 


Figure  4.  Mass  density  (kg/m3)  contours,  MHRDR  2-d 
calculation,  time  3.4  psec;  first  solid  line  adjacent  to 
dotted  lines  represents  4.55X10'2  kg/m3,  other  solid  lines 
each  20%  higher,  dotted  20%  lower. 
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